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1. INTRODUCTION 

The Prototype Fast Breeder Reactor (PFBR) is a 500 MW capacity sodium cooled fast breeder 
reactor under commissioning in Kalpakkam, India. The steam generator (SG) of PFBR is a counter flow shell 
and tube heat exchanger with sodium on shell side and water/steam on the tube side. The SG contains 547 
tubes of 17.2 mm OD, 25 m length and thickness of 2.3 mm. The wall thickness of the tube is the only barrier 
between the liquid sodium and the water in the SG. Periodic in-service inspection of the SG tubes is 
mandated to mitigate the failure of the SG tubes. The in-service inspection of the SG tubes is carried out 
using the indigenously developed PFBR SG inspection system (PSGIS) device. The structural integrity of the 
SG tubes is evaluated using the remote field eddy current testing technique. The PSGIS device consists of 
several modules such as the device deployment module, tube locator module (TLM), and cable pusher 
module. The primary task of the inspection device is to orient the cable pusher module above the user defined 
SG tube for the deployment and retrieval of the eddy current probe at uniform speed along the tube length. 
The TLM of the PSGIS device is a two-arm serial robotic device. It has two rotational degrees of freedom 
which enable it to position the end-effector of the robotic arm above anyone of the 547 tubes as shown in 
Figure 1. The actual prototype of two-link robotic arm is shown in Figure 2. 

Forward kinematics determines the end-effector position corresponding to joint angle. Figure 3 
represents the schematic of two-link robotic arm. The position of two-axis robotic arm (x, y) is represented 
by (1) and (2). 


x = 1, cos(@,) + L cos(@, + 82) (1) 
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y = L sin(84) + l, sin(@, + 62) (2) 


where x and y are end-effector position coordinates, 1; is the length of link1, lz is the length of link2, and 6; 
and 02 are joint angle of link! and link2, respectively. 
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Figure 1. SG tube sheet and steam generator 
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Figure 2. Prototype of two-link robotic arm Figure 3. Schematic of two-link robotic arm [1] 


TLM uses harmonic drive FHA-C miniseries (gear ratio=50) servo actuators with incremental line 
driver encoder (2000 ppr). Table 1 shows the two-link robotic arm geometric and mass parameter. Inverse 
kinematics establishes the relation between the joint angle and end-effector position. There are many 


methods to find out joint angles such as analytical, geometric, and iterative techniques. The inverse kinematic 
algorithm of two-link robotic arm is shown in Figure 4. 


Table 1. Two-link robotic arm geometric parameter 


Parameters Value Parameters Value 
Mass of the Elbow arm motor(ms3) 0.5 kg Elbow Arm Moment of Inertia about CG (Iz2) 3.675*10° kgm? 
Mass of the Elbow arm (m2) 1.0kg Length of shoulder arm (lı) 0.21 m 
Mass of the Shoulder Arm (m,) 1.0kg Length of Elbow arm (l2) 0.21 m 
Shoulder Arm Moment of Inertia about CGdz1) 3.675*10° kgm? Mass of the point load at end-effector end(m4) 1kg 
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Figure 4. Inverse kinematics algorithm of SG-ISI 


In TLM, damped least square inverse kinematics is used to find out joint velocity as given by (3), 
6 = FY + PAV) (3) 


where A is damping factor, V is end-effector velocity, g- joint velocity of two-link robotic arm, and J 

represents the Jacobian of the two-link robotic arm which is given by (4). 
=l sin(6,) = lL sin(6, + 02) =l; sin(6, + 02) 

J= (4) 

lL cos(81) + l, cos(6, +82) Ll, cos(@, + 02) 


The positioning accuracy required for the TLM is of the order of 0.1 mm. Moreover, in a serial 
manipulator like TLM, achieving accuracy of the order of 0.1 mm is a challenging task. The positioning 
accuracy is affected by residual vibration of the robotic arm. Hence, selection of smooth motion profile plays 
a vital role in achieving the 0.1 mm accuracy. The smooth motion profile is generated in many ways by 
different functions and different methods. 

Motion planning in robotics is a growing research field in the industry. Robotic motion is classified 
into two types, namely point-to-point and motion through sequence of points. Motion planning is done in 
joint space/end-effector coordinates. The end-effector coordinates, and joint coordinates are related by 
nonlinear functions, so planning of joint motion profile plays an important role in end-effector motion 
planning. In literature, several methods are proposed for the generation of smooth trajectory namely 
polynomial, spline, dynamic movement primitives, digital convolution, and finite impulse response filter- 
based methods. It is seen from the literature that Ezair et al. [2] and Ķirecci and Gilmartin [3] investigated 
motion planning from a computational point of view. The methods favorable for reducing the computation 
load are i) root multiplicity method, ii) recursive method, and iii) convolution method. 
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Ezair et al. [2] enumerated a recursive method to reduce higher-order polynomial to lower-order 
polynomial. It can produce a smooth trajectory using arbitrary initial and final conditions. Kirecci and 
Gilmartin [3] synthesized the joint angle motion planning using arbitrary power polynomial. The arbitrary 
power polynomial produces lower peak velocity than the other motion profile. Moreover, Boryga and Grabos 
[4] proposed root multiplicity method to form higher order polynomial. The higher-order polynomial 
produces oscillation. In addition, previous studies showed that higher-order polynomials introduce 
computational complexity. To avoid this, a piecewise polynomial (spline) is used to interpolate the joint 
angle [5]. Even though piecewise polynomial (spline) is used to join the number of via points for smooth 
motion, the acceleration of cubic spline is not smooth. Furthermore, truncated trigonometric series is used to 
interpolate the joint angle [6]. Moreover, motion planning is discussed from the point of minimum jerk [7]- 
[9], minimum time [10]—[12], minimum energy [13], [14], and minimum torque [15]. 

Jerk-limited motion is generated by polynomial, trigonometric and sigmoid function. M. S. Mujtaba 
[16] stated that bang-bang (quadratic polynomial) motion produces infinite jerk which produces lots of 
vibration which leads to tracking error. In addition, jerk-limited motion is generated by S-curve profile. S- 
curve motion profile has seven phases. Depending on the motion distance, initial and final boundary 
conditions, S-curve may not have constant velocity and constant acceleration phases. Nguyen et al. [17] 
proposed a recursive method to form S-curve motion profile using polynomial and trigonometric. They stated 
that position error decreases when higher-order polynomial S-curve is used. Further, an increase of 
polynomial degree increases the number of segments by 2"-1. The jerk-limited motion profile increases the 
motion time than cubic polynomial. 

At the beginning of the 19" century, only first-order continuity (velocity) is achieved in motion 
planning. At the end of the 19" century, second-order continuity (acceleration) is achieved. In the 20" 
century, third-order continuity (jerk) has been a main concern in the motion profile. Moreover, jerk-bounded 
motion produces symmetric and asymmetric motion profiles due to the bounded value of jerk in acceleration 
and deceleration phases. The design of acceleration and deceleration profile plays a vital role in improving 
precision and accuracy. Moreover, duration of acceleration and deceleration play a significant role in 
reducing the settling time and residual vibration. The settling time is reduced by asymmetric motion and 
lower magnitude of jerk. The conventional method cannot produce an asymmetric velocity profile. Different 
magnitudes of jerk in acceleration and deceleration can produce an asymmetric profile. Li [18] proposed an 
asymmetric motion profile using a sine jerk profile which reduces residual vibration. The vibration produced 
by the sinusoidal jerk profile is lesser than the trapezoidal and jerk-limited s-curve motion profile. Rew and 
Kim [19] enumerated closed-form solution to asymmetric motion using a seven-segment approach to reduce 
settling time and residual vibration. Amthor et al. [20] proposed an asymmetric motion profile using fourth- 
order polynomial. Perumal and Jawahar [21] proposed a synchronized trigonometric S-curve for jerk 
reduction. Moreover, sigmoid (piecewise exponential) is also used for smooth motion [22], [23]. 

The motion profile with higher-order continuity increases computational complexity but motion 
time is more complex, and it needs precise hypothesis. Snap bounded trajectory planning reduces the 
amplitudes of vibration than the jerk-limited motion profile. Rocha et al. [24] proposed an embedded system- 
based snap constrained trajectory planning method for 3D motion systems. They compared snap constrained 
trajectory planning with seven segments and trapezoidal acceleration motion profile. 

Even though piecewise polynomial (algebraic spline) is used to plan the joint motion, it is not 
suitable for operations such as obstacle avoidance. Moreover, drawback of algebraic spline is that 
computational complexity increases due to coupled nature of unknowns. To reduce the coupling, Simon et al. 
[25] introduced a trigonometric spline to interpolate data points. Further, Visioli [26] compared the algebraic 
and trigonometric splines to reduce the overshot of the motion profile. Nnaji and Asano [27] investigated five 
joint trajectory profiles such as bang-bang, polynomial, exponential, cosine, and sine ramp profile. They 
investigated joint motion based on total work done, individual peak joint torque, and individual peak joint 
power. They stated that the evaluation process is independent of joint profile or specific robot. They stated 
that exponential trajectory was the best path out of the five trajectories. Polynomial trajectory is robot 
specific but exponential trajectory is not robot specific. The exponential profile is more suitable for large 
distance motion than the shorter motion. The acceleration produced by the exponential motion profile at the 
starting and ending is finite which results in tracking error. 

On contrary to off-line programming, online motion planning needs fewer data points. Moreover, 
online motion planning plays a vital role in avoiding obstacles. The finite impulse response filter is used to 
generate multi-segment motion profile in online motion planning [28]. Chand and Doty [29] proposed a cubic 
spline online joint trajectory. The average path error produced by the cubic spline is lower than the linear 
spline joint trajectory. Bazaz and Tondu [30] proposed a three-cubic method for motion generation between 
two successive points. Sidobre and Desormeaux [31] proposed a smooth cubic polynomial for human-robot 
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interactions. In robotics motion control, shorter movements possess complex behavior. They analyzed the 
influence of velocity on the shorter motion. 

From the literature, it was found that no generalized study was carried out for smooth motion profile 
with kinematic constraints. In this paper, generalized study of bang-bang, jerk-limited s-curve, fifth order 
polynomial (quintic polynomial), and cycloid motion profile is presented. The paper is organized as follows. 
Section 1 gives a brief review of forward, inverse kinematics of two link robotic arm and state of the art in 
motion planning. Section 2 gives a brief overview of the smooth motion profiles. Section 3 gives analytical 
methods to compute peak velocity, acceleration, and jerk. Section 4 discusses case study. Section 5 
summarizes the results of this paper. 


2. SELECTION OF SMOOTH MOTION PROFILE FOR STEAM GENERATOR INSPECTION 

Smooth motion is obtained in many ways. Smooth motion is generated by with or without 
considering kinematic and dynamic constraints. The kinematic constraints are joint actuator velocity, 
acceleration, and jerk. The dynamic constraints are joint actuator torque and torque rate. Further, minimum 
time is handled by kinematic and dynamic constraints, respectively. The bang-bang motion profile achieves 
time-optimal motion, but it results in jerky motion due to impulse load (infinite jerk) at finite points as shown 
in Figure 5(a), which causes tracking error and vibration. However, the jerky motion is smoothened by linear 
acceleration and deceleration (trapezoidal acceleration and deceleration) which increases motion time as 
shown in Figure 5(b). Even though jerk-limited motion produces finite jerks, it does not have continuous 
jerks as shown in Figure 5(c). This can be further improved by including higher-order continuity namely 
derivative of jerk (snap), derivative of snap (jounce or crackle), and which makes the jerk continuous. The 
motion profile with higher-order continuity results in computational complexity which needs proper 
hypothesis for selecting motion time. The higher-order polynomial motion produces smaller jerks. The trade- 
off has to be made such that motion profile satisfies computational complexity and bounded jerk. In multi- 
segment motion planning, the movement from one segment to another segment needs smooth transition 
(preferably c? continuity) to avoid impulse load (impulse force or torque). The piecewise polynomial is used 
to join the segment smoothly. The formulation of various motion profiles is given in the methodology 
section. 
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Figure 5. Position, velocity, acceleration, and jerk of various motion profiles: (a) bang-bang, (b) cubic 
polynomial, (c) jerk-limited s curve, (d) fifth-order polynomial, and (e) cycloid motion 


3. METHODOLOGY 
In this section, the end-effector peak velocity, peak acceleration and peak jerk of various motion 
profiles are derived and compared for the given motion time (t) and total displacement (As). 
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3.1. Quintic polynomial 

Polynomial is used for motion planning of joint angle due to less computation. Even though linear 
polynomials have position continuity, they do not have velocity, acceleration continuity. On the other hand, 
cubic polynomial has position and velocity continuity. However, the acceleration produces discontinuity at 
the starting and ending in cubic polynomial. The fifth-order polynomial (Figure 5(d)) has position, velocity, 
and acceleration continuity. Moreover, jerk produced by the cubic polynomial is constant, but fifth-order 
polynomial produces non-linear jerk. The position (s), velocity (v), acceleration (a), jerk (j), and snap of fifth 
order polynomial is given by (5), 


s(t) = ao + at + azt? + azt? + a,t* + ast* 
v(t) = a, + 2a,t + 3a3t? + 4a,t? + 5ast* 
a(t) = 2a, + 6a3t + 12a,t? + 20a;t? (5) 
j(t) = 6a3 + 24a,t + 60ast? 
j(t) = 24a, + 120ast 
where ao,a1,a2,a3,a4,a5 are unknowns and t-time. 
The fifth-order polynomial has six unknowns. The six unknowns are determined by applying the 


initial and final boundary conditions. The initial and final boundary conditions are given by (6). 


s(t=0)=s; vit=0)=0 at=0)=0 


s(t=t,)=s, v(t=t)=0 a(t=t,)=0 (6) 


The position at time t is given by (7) 


s(t) = (G t? — (=) ay (=) | (sp — si) (7) 


The jerk of fifth-order polynomial is given by (8) 


o- -eea » 


The maximum velocity occurs when acceleration is zero (t=t/2).The maximum velocity (Vm) of fifth-order 
polynomial is given by (9a) and (9b) 


m= È) ov 


whereAs = sp — Si 
Vm = 0.4789(As? jn) 2) (9b) 


where vmis the maximum velocity, tr is the motion time, As is the total distance, and jm is the maximum jerk. 
The maximum value of jerk occurs at t=0 and t=t;. The given distance (As) and velocity (vm), the maximum 
jerk Gm) is calculated by (10) 


3 


am (10) 


Jm = (As)2(0.4789)3 


where jm is the maximum jerk, Vm is the maximum velocity, and As is the total distance 
The optimal time for fifth-order polynomial is given by (11) 


as\ 1/3) 
ty = 3.9149 (=) (11) 
The maximum acceleration (am) of fifth-order polynomial is given by (12) 
am = 0.3767(Asj2,)°/® (12) 
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3.2. Jerk-limited S-curve (without constant acceleration and constant velocity phase) 

The cubic polynomial produces linear acceleration, but this introduces finite jerks at the starting and 
ending. Even though the jerk-limited motion profile has some finite value at the starting and ending in jerk, it 
has zero acceleration at endpoints. Depending on the distance and kinematic constraints (maximum 
acceleration and maximum velocity), some phases (constant acceleration or constant velocity) may not exist. 
For ease of computational complexity, only constant jerk phase is considered. The constant acceleration and 
constant velocity phase are not considered. It can be seen that the minimum number of segments for 
achieving minimum time is 3 for satisfying only jerk constraints without utilizing the maximum acceleration 
and velocity. 

The jerk (j) at time t is given by (13), 


jm 0<t<tj 
jm 3t < t< 4t; 


where tj-time duration for jerk phase, jm- maximum jerk. The position of jerk-limited motion profile is given 
by (14a). 


jm (=) 0<t<t 
O5] M-A an) uessa gai 
MORTODATO ROEE <0 44 
The velocity of jerk-limited motion profile is given by (14b). 


. (t? 
jm (=) 0<t< tj 
; . (t? . (F 
v(t) = Jmtjt —Jm (=) + Jm (+) 0<t< 3t (14b) 


a ft? ; . (ï 
jm (=) — jmtjt + jm () 3t; < t < 4tj 


The acceleration of jerk-limited motion profile is given by (14c). 


jmt 0<t<t 
jm 3t; <t <4t, 


The time duration for jerk phase along x, y-direction (tx, ty) is given by (15). 


r-T” = fy (15) 


From the maximum velocity, the maximum jerk is calculated. From the maximum jerk (jm) and distance, 
motion time is calculated. The time for the jerk phase is given by (16a). 


t; = 0.7937 a (16a) 
The total time for jerk-limited motion profile is given by (16c). 

tj = 4t; 

t; = 3.1748 (ayn (16b) 

to = max(ty, ty) (16c) 
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3.3. Cycloid motion 

The fifth-order polynomial and jerk-limited motion are formed by polynomial, whereas cycloid 
motion (Figure 5(e)) is constructed by point p on the circle rolls along a straight line as shown in Figure 6. 
The cycloid motion cannot be used for motion planning of multiple segments, due to finite jerk in the 
intermediate points, start and endpoints. The cycloid motion has linear polynomial and trigonometric 
function. 


(sf-si) 


Figure 6. Generation of cycloid motion 


The position of cycloid motion profile as shown in Figure 6 is given by (17) . 


s(t) = As {(é) — (=) sin (=) + 5; (17) 


The velocity, acceleration and jerk are given by(18-20). 


v(t) = (=) {1 — cos (=)} (18) 
a(t) = (F) {sin (2) (19) 
w- (2) om) an 


where As is the total displacement, tr is the total motion time. The acceleration profile of cycloid motion is 
smoother than the jerk-limited motion profile. The peak acceleration of cycloid motion is larger than the 
fifth-order polynomial. 


3.4. Harmonic jerk S-curve 
The harmonic jerk S-curve is given by (21). 


im (1 = cos (##)) 0<t<tj 

—jm (: — COS e) ti <ts2t; 
—jm (: — cos =) 2t; <t<3tj 
da (: — COS e=) 3t; <t < 4t; 


The acceleration of harmonic jerk S-curve is given by (22). 


j&© = (21) 
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0<t<tj 


—0.5 * jm 4t — (=) +054 jm*t) t)<t< 2t 
gi 
as J 
a(t) i : (=) (22) 
sın T 
—0.5 * jm 4 t = 2t; <t < 3t 
(5) 
sin( 222) 
0.5 * jm 4t — Lf —0.5%jm*t, 3 <t<4t; 


.\2 
+ jm * 0.5% tjt + jm *0.25*t? + jm*05*(L)  t)<t<2t 


2m 


v(t) = (23) 
2 cos( 2 t; 2 
—0.5 * jm E Ee + jm * 0.5 * tÊ + jm * 0.5 * (=) 2t; <t <3t; 
(5) 
amt 
5 2 cos(#) : 2 3 tj 2 2 
0.5 * jm at on + jm * 0.25 * EF = jm * 0.5 * (+) —jm * 0.5 * tit 3tj <t<4t; 
E) 
The position of harmonic jerk S-curve is given by (24). 
oi D) sas 0<t<t 

oil , a bint 05+ (H ein +025+Hf q+ 0250507 tine (2) EET geese 

st) = 7 (24) 


3 sin( 224) ; ee 
0.5 sf fs) H jm * 0.5 * tPt + (£) * t(=) + jm *0.5 xt? 2tj < t< 3t 
5) 
y 2nt 
sin( 4) 


E l-n * 0.25 * tjt? + jm * 0.25 * tft (*) *t* (2) bind ee i g , (6) xt (2) ee 
‘ty, 


4. RESULTS AND DISCUSSION 

The position of two link robotic arm from start point to end point is shown in Table 2. It is assumed 
that end-effector velocity along x-direction is 20 mm/sec. Bang-bang motion produces finite acceleration at 
starting and ending, which makes the motion jerkiness. On the contrary, cubic polynomial produces finite 
jerk as shown in Table 3. The jerk-limited motion profile gives maximum velocity than other motion profiles 
for the given jerk (j) and motion time (ts). However, quintic polynomial produces the lowest velocity and 
acceleration. The maximum velocity and acceleration produced by S-curve is lower than the other motion 
profile. The comparison of different motion parameters is given in Table 3. 

The motion of x-coordinate of end-effector is shown in Figure 7. Figure 8 represents the motion of 
y-coordinate. The total time for both coordinates is the same. For fixed time, the maximum jerk magnitude 
along x-coordinate differs from the y-coordinate due to variation in displacement in the x and y coordinates. 
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Table 2. End-effector displacement at the start and end point 


End-effector displacement Initial point (mm) End point (mm) Kinematics constraints 
x 209.3 225.4 eee 
¥ 250.9 223 1 End-effector velocity along x-direction= 20 (mm/sec) 
Table 3. Comparison of different motion parameters for the given time 
Type of Interpolation Total distance (mm) J (mm/s?) a (mm/s?) v (mm/sec) t(s) 
Bang-bang motion As co am af Gm As ty 
Cubic polynomial As Timi: 1.15 (Jinas 0.66 As Zjm1) tf 


Jerk-limited motion 
Quintic polynomial 
Cycloid motion 


Harmonic jerk S-curve 


As Tho 0.79 


) 066( 
) 065( 

As jms 0.38 (; As 0.48 ({/As7jina) ty 
(ese)  os9 


As Fa 0.54 


As Jms 0.5 (; fass) 0.5 (/AsZjms) ty 


End effector x-coordinate displacement(mm) 


— Jerk limited S-curve 
~--Quintic polynomial 
~--Cycloid 

Harmonic jerk limited S-urve' 
* Cubic polynomial 


L 


0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 
Time(sec) 


Figure 7. End-effector x-coordinate with respect to time 
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Figure 8. End-effector y-coordinate with respect to time 
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It can be seen that only time is synchronized, while velocity is not synchronized as shown in Figures 9 
and 10. The maximum velocity of x-coordinate differs from y-coordinate velocity due to different 
displacement along x and y coordinates. The peak velocity of cubic polynomial is lower than the other 
motion profile. It can be seen that cycloid motion produces higher peak velocity than quintic motion profile. 
However, quintic polynomial produces the lowest peak velocity with smoothness. 
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Figure 9. End-effector x-coordinate velocity with respect to time 
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Figure 10. End-effector y-coordinate velocity with respect to time 


The acceleration produced by the jerk-limited motion is not smooth. Moreover, acceleration is linear 
in jerk-limited motion whereas it is cubic in quintic polynomial. In contrast to cubic polynomial, acceleration 
of cycloid motion is zero at the starting and ending but the jerk produced by the cycloid motion is finite at the 
starting and ending which leads to trajectory error. The acceleration produced by the jerk-limited S-curve is 
linear as shown in Figures 11 and 12. 

It can be seen that jerk produced by the S-curve is a step function. In quintic and cycloid motion 
profile jerk is continuous except at the starting and ending (finite jerk) unlike step function in jerk-limited S- 
curve as shown in Figure 13. The jerk magnitude of cycloid motion is higher than the jerk-limited S-curve as 
shown in Table 4. The maximum jerk value in the x and y direction is different due to synchronization of 
time as shown in Figures 13 and 14. For time synchronized motion, jerk in the x-direction is lesser than the y- 
direction. It was found that harmonic jerk S-curve produces maximum jerk for the given motion time as 
shown in Table 4. 

Figure 15 shows the end-effector profile of a two-link robotic arm. The end-effector coordinates are 
interpolated by fifth-order polynomial. The position of joint angle for the initial and final point is calculated 
by inverse kinematics algorithm. 
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Figure 11. End-effector x-coordinate acceleration with respect to time 
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Figure 12. End-effector y-coordinate acceleration with respect to time 
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Figure 13. End-effector x-coordinate jerk with respect to time 
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Table 4. Comparison of jerk magnitude end-effector x-coordinate for various motion planning 


Motion model Maximum jerk value (mm/s?) Maximum acceleration (mm/s?) Maximum velocity (mm/s) 
Cubic polynomial 46.3 -46.3 15 

Jerk-limited S-curve 123.43 49.68 20 

Harmonic jerk S-curve 246.90 49.68 20 

Quintic polynomial 231.43 35.85 18.75 

Cycloid motion 152.27 39.02 20 
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Figure 14. End-effector x-coordinate jerk with respect to time 


300 T T T T T T 


a O O O oO O O Q O 
start point 
250 © Q 


2000 


Y(mm) 


150F 
O 


100 F 


50 f 


0 50 100 150 200 250 300 
X(mm) 


Figure 15. Two-link robotic arm end-effector profile for fifth-order polynomial 


5. CONCLUSION 

In this paper, the motion analysis of two-link robotic arm for steam generator inspection is discussed 
from kinematics aspects without considering dynamics. The theoretical formulation of fifth-order, jerk- 
limited and cycloid motion profile is presented, and optimal time calculation has been presented for the 
velocity constraints. Finally, a case study has been presented. The conclusions are summarized as follows: i) 
the end effector trajectory is planned by using fifth-order polynomial, cycloid motion, harmonic jerk S-curve 
and jerk-limited S-curve, ii) the time-optimality of motion profile considering velocity and acceleration 
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constraints is discussed, and iii) the synchronization of end-effector coordinates motion along the x and y 
direction is discussed. In future, this study can be extended to motion planning with obstacles in between the 
start point and end point. Moreover, comparison of various motion profiles with dynamic constraints will be 
studied. 


ACKNOWLEDGEMENTS 

Support and encouragement provided by Shri S. Raghupathy, Director, RDTG, IGCAR, and Dr. 
A.K. Bhaduri, Director, IGCAR are thankfully acknowledged. Perumalsamy G would like to thank 
Department of Atomic Energy fellowship for a perspective research grant for Ph.D. 


REFERENCES 


[1] G. Perumalsamy, P. Visweswaran, J. Jose, S. Joseph Winston, and S. Murugan, “Quintic interpolation joint trajectory for the path 
planning of a serial two-axis robotic arm for PFBR steam generator inspection,” in Machines, Mechanism and Robotics, 2019, pp. 
637—648. 

[2] B. Ezair, T. Tassa, and Z. Shiller, “Planning high order trajectories with general initial and final conditions and asymmetric 
bounds,” The International Journal of Robotics Research, vol. 33, no. 6, pp. 898-916, May 2014, doi: 
10.1177/0278364913517148. 

[3] A. Kirecci and M. J. Gilmartin, “Improved trajectory planning using arbitrary power polynomials,” Proceedings of the Institution 
of Mechanical Engineers, Part I: Journal of Systems and Control Engineering, vol. 208, no. 1, pp. 3-13, Feb. 1994, doi: 
10.1243/PIME_PROC_1994_208_301_02. 

[4] M. Boryga and A. Graboś, “Planning of manipulator motion trajectory with higher-degree polynomials use,” Mechanism and 
Machine Theory, vol. 44, no. 7, pp. 1400-1419, Jul. 2009, doi: 10.1016/j.mechmachtheory.2008.11.003. 

[5] Y. Fang, J. Hu, W. Liu, Q. Shao, J. Qi, and Y. Peng, “Smooth and time-optimal S-curve trajectory planning for automated robots 
and machines,” Mechanism and Machine Theory, vol. 137, pp. 127—153, Jul. 2019, doi: 10.1016/j.mechmachtheory.2019.03.019. 

[6] M. L. Hornick and B. Ravani, “Computer-aided off-line planning and programming of robot motion,” The International Journal 
of Robotics Research, vol. 4, no. 4, pp. 18-31, Jan. 1986, doi: 10.1177/027836498600400402. 

[7] Y. Fang, J. Qi, J. Hu, W. Wang, and Y. Peng, “An approach for jerk-continuous trajectory generation of robotic manipulators with 
kinematical constraints,” Mechanism and Machine Theory, vol. 153, Nov. 2020, doi: 10.1016/j.mechmachtheory.2020.103957. 

[8] A. Piazzi and A. Visioli, “Global minimum-jerk trajectory planning of robot manipulators,’ IEEE Transactions on Industrial 
Electronics, vol. 47, no. 1, pp. 140-149, 2000, doi: 10.1109/41.824136. 

[9] H.-L Lin, “A fast and unified method to find a minimum-jerk robot joint trajectory using particle swarm optimization,” Journal of 
Intelligent & Robotic Systems, vol. 75, no. 3—4, pp. 379-392, Sep. 2014, doi: 10.1007/s10846-013-9982-8. 

[10] F. J. Abu-Dakka, I. F. Assad, R. M. Alkhdour, and M. Abderahim, “Statistical evaluation of an evolutionary algorithm for 
minimum time trajectory planning problem for industrial robots,” The International Journal of Advanced Manufacturing 
Technology, vol. 89, no. 1—4, pp. 389-406, Mar. 2017, doi: 10.1007/s00170-016-9050-1. 

[11] D. Constantinescu and E. A. Croft, “Smooth and time-optimal trajectory planning for industrial manipulators along specified 
paths,” Journal of Robotic Systems, vol. 17, no. 5, pp. 233-249, May 2000, doi: 10.1002/(SICI)1097- 
4563(200005)17:5<233:AID-ROB1>3.0.CO;2-Y. 

[12] A. Piazzi and A. Visioli, “Global minimum-time trajectory planning of mechanical manipulators using interval analysis,” 
International Journal of Control, vol. 71, no. 4, pp. 631-652, Jan. 1998, doi: 10.1080/002071798221713. 

[13] G. Berselli, F. Balugani, M. Pellicciari, and M. Gadaleta, “Energy-optimal motions for servo-systems: A comparison of spline 
interpolants and performance indexes using a CAD-based approach,” Robotics and Computer-Integrated Manufacturing, vol. 40, 
pp. 55-65, Aug. 2016, doi: 10.1016/j.rcim.2016.01.003. 

[14] P. Boscariol and D. Richiedei, “Energy-efficient design of multipoint trajectories for Cartesian robots,” The International Journal 
of Advanced Manufacturing Technology, vol. 102, no. 5-8, pp. 1853—1870, Jun. 2019, doi: 10.1007/s00170-018-03234-4. 

[15] T. Chettibi, H. E. Lehtihet, M. Haddad, and S. Hanchi, “Minimum cost trajectory planning for industrial robots,” European 
Journal of Mechanics - A/Solids, vol. 23, no. 4, pp. 703-715, Jul. 2004, doi: 10.1016/j.euromechsol.2004.02.006. 

[16] M. S. Mujtaba, “Discussion of trajectory calculation methods,” Exploratory Study of Computer Integrated Assembly Systems, 
California, 1977. 

[17] K. D. Nguyen, T.-C. Ng, and I.-M. Chen, “On algorithms for planning S-curve motion profiles,” International Journal of 
Advanced Robotic Systems, vol. 5, no. 1, pp. 99-106, Mar. 2008, doi: 10.5772/5652. 

[18] H. Li, “A jerk-constrained asymmetric motion profile for high-speed motion stages to reduce residual vibration,” International 
Journal of Computer Applications in Technology, vol. 53, no. 2, pp. 149-156, 2016, doi: 10.1504/IJCAT.2016.074453. 

[19] Keun-Ho Rew and Kyung-Soo Kim, “A closed-form solution to asymmetric motion profile allowing acceleration manipulation,” 
IEEE Transactions on Industrial Electronics, vol. 57, no. 7, pp. 2499-2506, Jul. 2010, doi: 10.1109/TIE.2009.2036032. 

[20] A. Amthor, J. Werner, A. Lorenz, S. Zschaeck, and C. Ament, “Asymmetric motion profile planning for nanopositioning and 
nanomeasuring machines,” Proceedings of the Institution of Mechanical Engineers, Part I: Journal of Systems and Control 
Engineering, vol. 224, no. 1, pp. 79-92, Feb. 2010, doi: 10.1243/095965 18JSCE826. 

[21] S. Perumaal and N. Jawahar, “Synchronized trigonometric s-curve trajectory for jerk-bounded time-optimal pick and place 
operation,” International Journal of Robotics and Automation, vol. 27, no. 4, pp. 385-395, 2012, doi: 
10.2316/Journal.206.2012.4.206-3780. 

[22] C. DiMarco, J. C. Ziegert, and C. Vermillion, “Exponential and sigmoid-interpolated machining trajectories,’ Journal of 
Manufacturing Systems, vol. 37, pp. 535-541, Oct. 2015, doi: 10.1016/j.jmsy.2015.04.007. 

[23] Z. Rymansaib, P. Iravani, and M. N. Sahinkaya, “Exponential trajectory generation for point to point motions,” in 20/3 
IEEE/ASME International Conference on Advanced Intelligent Mechatronics, Jul. 2013, pp. 906-911, doi: 
10.1109/AIM.2013.6584209. 

[24] P. A. S. Da Rocha, W. D. De Oliveira, and M. E. De Lima Tostes, “An embedded system-based snap constrained trajectory 
planning method for 3D motion systems,” IEEE Access, vol. 7, pp. 125188—125204, 2019, doi: 10.1109/ACCESS.2019.29391 16. 


IAES Int J Rob & Autom, Vol. 11, No. 3, September 2022: 181-195 


IAES Int J Rob & Autom ISSN: 2722-2586 Oo 195 


[25] 
[26] 
[27] 
[28] 
[29] 


[30] 


D. Simon and C. Isik, “Optimal trigonometric robot joint trajectories,” Robotica, vol. 9, no. 4, pp. 379-386, Dec. 1991, doi: 
10.1017/S0263574700000552. 

A. Visioli, “Trajectory planning of robot manipulators by using algebraic and trigonometric splines,” Robotica, vol. 18, no. 6, pp. 
611-631, Nov. 2000, doi: 10.1017/S0263574700002721. 

B. O. Nnaji and D. K. Asano, “Evaluation of trajectories for different classes of robots,” Robotics and Computer-Integrated 
Manufacturing, vol. 6, no. 1, pp. 25-35, Jan. 1989, doi: 10.1016/0736-5845(89)90082-3. 

L. Biagiotti and C. Melchiorri, “Trajectory generation via FIR filters: A procedure for time-optimization under kinematic and 
frequency constraints,” Control Engineering Practice, vol. 87, pp. 43-58, Jun. 2019, doi: 10.1016/j.conengprac.2019.03.017. 

S. Chand and K. L. Doty, “On-line polynomial trajectories for robot manipulators,” The International Journal of Robotics 
Research, vol. 4, no. 2, pp. 38—48, Jun. 1985, doi: 10.1177/027836498500400204. 

S. A. Bazaz and B. Tondu, “Minimum time on-line joint trajectory generator based on low order spline method for industrial 
manipulators,” Robotics and Autonomous Systems, vol. 29, no. 4, pp. 257—268, Dec. 1999, doi: 10.1016/S0921-8890(99)00058-5. 


BIOGRAPHIES OF AUTHORS 


G Perumalsamy © ki P obtained his M.Des. in Mechanical Systems in IITDM, 
Kancheepuram. He is currently pursuing his doctoral degree in Engineering Science at HBNI, 
Mumbai, India. His thesis focuses on smooth trajectory profile generation for steam generator 
inspection in nuclear reactor. E-mail: perumalsamy2012mdes @ gmail.com. 


P Visweswaran © EJ P' has been working as Scientific Officer at Indira Gandhi Centre 
for Atomic Research (IGCAR), Kalpakkam, India. His areas of interest are remote handling 
and robotic development. E-mail: visweswaran @igcar.gov.in. 


a 


a i 
Deepak Kumar Ok P has been working at IGCAR, Kalpakkam, India since 2017. His 
areas of interest are remote handling and robotic development. E-mail: deepakk @igcar.gov.in. 
A f 


S. Joseph Winston EJ P has been working as Scientific Officer and The Head Remote 
Handling and Irradiation Experiments Division of Indira Gandhi Centre for Atomic Research 
(IGCAR), Kalpakkam, India. His areas of interest are in the development of remote inspection 
systems. E-mail: winston @igcar.gov.in. 


=— S. Murugan © ki P' was with Reactor Design and Technology Group of IGCAR. His 
areas of interest are irradiation experiments, numerical modeling in welding, residual stress 
analysis, precision machining and fabrication, robotics and remote tooling, in-service 
inspection, and fast reactor engineering development. E-mail: sankmurugan1 @ gmail.com. 


Selection of smooth motion profile for a tube locator module of an inspection device (G. Perumalsamy) 


